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Abstract. A family of explicit Lyapunov functions for positive recurrent Mar- 
kovian Jackson networks is constructed. With this result we obtain explicit 
estimates of the tail distribution of the first time when the process returns 
to large compact sets and some explicit estimates of the essential radius of 
the process. The essential spectral radius of the process provides the best 
geometric convergence rate to equilibrium that one can get by changing the 
transitions of the process in a finite set. 



1. Introduction 

Before formulating our results we recall the definition and some well known re- 
sults concerning classical Jackson networks, see [7] for example. For a Jackson 
network with d queues, the arrivals at the i-th queue are Poisson with parameter \ 
and the services delivered by the server are exponentially distributed with parame- 
ters [li. All the Poisson processes and the services are assumed to be independent. 
The routing matrix is denoted P — (pij] i,j = 1, . . . , d), pij is the probability that 
a customer goes to the j'-th queue when he has finished his service at queue i. The 
residual quantity 

d 

PiO = 1 - ^2 Pij 
3=1 

is the probability that this customer leaves definitively the network. Without any 
restriction of generality we can assume that pa — for all i € { 1, . . . , d} . 

Denote by Zi(t) the length of the queue i at time t. Then the process Z(t) = 
(Zi, (t), . . . , Zd{t)) is a continuous time Markov process on generated by 

CM = E *)(/(*) - f(y)), y G <, 

with q(y, z) = q(z — y) such that 

Xi, ify = e l , ie {1, . ..,d}, 
ViPiO, if V = -e 4 , i G {1, • ■ • ,d}, 
(HPij, iiy = e J - e\ i,j E {1, . . .,d}, 



(i-i) q(y) 



0, otherwise, 



where e l denotes the ith unit vector, = if j ^ i and e\ = 1. It is convenient to 
put poo = 1 and Poi = for i ^ 0, the matrix (jpij] i, j — 0, ... ,d) is then stochastic. 
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(n) 

We denote by p\^ the n-time transition probabilities of a Markov chain with d + 1 
states associated to the stochastic matrix (jpif, i,j — 0, . . . ,d). 

Assumption (A). We suppose that the matrix (q(x — y); x,y G Z d ) is irre- 
ducible. 

This assumption is equivalent to the following conditions 

(A\) Every customer leaves the network with probability 1, i.e. for any i G 
{1, . . . , d} there exists neN, such that p^ > 0. This condition is satisfied 
if and only if the spectral radius of the matrix (jpij ; i, j = 1 , . . . , d) is strictly 
less than unity. 

(A 2 ) for any i = l,...,d, there exist n G N and j G {1, . . . , d} such that XjP^ > 
0. 

Under the assumption (A\), the system of traffic equations 

d 

(1.2) Vj = A 3 ' + 22^iPij, Q=l,..., d. 

i=i 

has a unique solution (z/j), and this solution satisfies Vi > for all i G {1, . . . , d}. 
The Markov process (Z(t)) is ergodic (positive recurrent) if and only if 

(1.3) Vi < fii for all % = 1, . . . , d, 

and the stationary probabilities (tt(x); x € Zi) are given by the product formulae 

d 

(1.4) w(x)=Y[(v i / fM ) Xi {l-v i / fM ), xeZ d + . 

»=i 

Assumption (B). We assume that the inequalities (|1.3p hold. 

Fayolle, Malyshev, Men'shikov and Sidorenko [3] proved that the rate of conver- 
gence to stationary distribution for ergodic Jackson networks is exponential. The 
proof of this result relies on the construction of a positive Lipschitz continuous 
function / : — > [0, +00 [ satisfying the inequality 

(1.5) Cf(x) < -e, yxeZl\E 

for some e > and some finite subset E C lA. Such a function / is often called a 
Lyapunov function for the Markov process (Z(t)). Using (|1.5j) one can easily show 
that for a > small enough, the function h(x) = exp(cr/(x)) satisfies the inequality 

(1.6) Ch{x) < -6h(x), VxeZ d + \E 

for some = 8(e) > 0. Usually, a function h : — > R + satisfying the inequality 
(|1.6p and such that 

(1.7) c E (f) d ^ inf f{x) > 

is also called a Lyapunov function for (Z(i)). To make a difference with a Lyapunov 
function satisfying the inequality (|1.5j) . we call such a function h a multiplicative 
Lyapunov function. For the hitting time r^; = inf{i > : Z(t) G -E}, the inequali- 
ties (HHJ) and (|T7| imply that 

(1.8) V x {t e >t)< —^—E x (h(Z(t)), T E >t) < —^exp(-6t)h(x), 
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for all x £ Z 1 ^ \ E. An explicit form for the multiplicative Lyapunov function h and 
the quantity 9 whould therefore imply explicit estimates for the tail distribution 
of the hitting time Te. Unfortunately, construction of an explicit multiplicative 
Lyapunov function satisfying (jl.6l) for a given finite set E C Zl with the best 
possible 9 is usually a very difficult problem. In [3] , the Lyapunov function / itself 
and the corresponding set E are both rather implicit. 

In the present paper we construct a class of explicit multiplicative Lyapunov 
functions h : — > [1, +oo[ with an explicit 

6h d = — \im sup Ch(x)/h(x) > 0. 

|a;|^oo 

For any such a function h and any < 9 < 9^, one could therefore identify the 
set E where (|1.6p holds and get an explicit estimate for the tail distribution of the 
hitting time te- 

Using the explicit form of the Lyapunov functions we obtain an explicit estimate 
for the essential spectral radius of the process (Z(t)). Recall that the spectral radius 
r* of the process (Z(t)) is defined as the infimum of all those r > for which 

r^P x (Z(t) =y)dt< +oo, Vx,y E Z d + . 

o 

When the process (Z(t)) is recurrent we obviously have r* = 1. The essential 
spectral radius r* of (Z(t)) is the infimum of all those r > for which there is a 
finite set E C Z+ such that 

j* oo 

/ r"*P x (Z(t) =y, t e > t) dt < +oo for all x,y eZ d + \E. 
Jo 

For the recurrent Markov process (Z(t)), the quantity r* is equal to the infimum 
of all those r > for which there is a finite set E c 1>\ such that 

(1.9) / r- l P x (t e > t) dt < +oo for all x£Z d + \E 



o 



(see for instance Proposition 3.6 of [B]). Remark that for those r > for which 
(|1.9[) holds, the function 

i rs d£ f \ XT ' r-'V^TE > t) dt, for xeZ d + \E, 
• V ' [0 for xeE, 

satisfies the inequalities (|1.6p and (|1.7p with a given E, 9 = — logr and 

poo / d \ _1 

c E (f)>J r- t e-^ x '+^ t dt > hnr + ^Xi+fii)] 

The last property of the essential spectral radius r* combined with the estimates 
(|1.8p shows therefore that the quantity 9* e = — log r * is equal to the supremum of 
all 9 > for which there exists a multiplicative Lyapunov function h : TL d + — > K + 
satisfying the inequalities (|1.6|) and (jTTTJ) for some finite subset _E C Zl. This is 
also the best 9 > one could expect to have in (|1.8p . 

The essential spectral radius is moreover related to the rate of convergence to 
equilibrium. To calculate the rate of convergence to equilibrium, one should identify 
the spectral gap of the transition operator, and except for some very particular 
processes, this is an extremely difficult problem. Explicit estimates of the rate of 
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convergence are therefore of interest. Malyshev and Spieksma [8j proved that for 
some general class of Markov chains, the quantity r* gives an accurate bound for 
that : this is the best geometric convergence rate one can get by changing the 
transitions of the process on finite subsets of states. By Perssons principle (see 
Liming Wu |10j ) , for symmetric Markov chains the quantity r* is related to the in- 
essential spectral radius of the corresponding Markov semi-group. For more details 
concerning the relationship between the quantity r* and the rate of convergence to 
equilibrium see Liming Wu |10j . 

In [6 , the quantity r * was represented in terms of the sample path large deviation 
rate function I[o.t](') of the scaled processes Z e {t) = eZ(t/e), t e [0, T]. Recall 
that the family of scaled Markov processes (Z e (t), t E [0,T]) satisfies the sample 
path large deviation principle (see [TJ [2j [5j 0] ) with a good rate function T^tjG)- 
Corollary 7.1 of the paper [6] proves that 

(1.10) logr e * = - inf I[oM 

where the infimum is taken over all absolutely continuous functions cf> : [0, 1] — > M.i_ 
with cj)(0) = (f)(1) and such that 4>(t) ^ for all < t < 1. For d < 2, the quantity 
r* was calculated explicitly : in this case, the infimum at the right hand side of 
(|1.10[) is achieved at some constant function 4>(t) = x = (x 1 , . . . ,x d ) e with 
x % > for some 1 < i < d and x^ = for j ^ i. For d — 1, Proposition 7.1 of [5] 
shows that 



logr-: = -(^l-^xlf 
and by Proposition 7.2 of 6 , for d = 2, 

(1.11) logr* = -(1 -pi2P2i)min{(VMi- V^i) 2 ,(VJ^2" V^f}- 

Unfortunately, for higher dimensions d > 3, the variational problem (|1.10p seems 
very difficult to resolve. In the present paper, using the explicit Lyapunov functions, 
we obtain explicit estimates for the essential spectral radius r* for an arbitrary 
dimension d. The quantity r* is calculated explicitly for several examples of Jackson 
networks. 



2. General results 
To formulate our results, we need to introduce some additional notation : 

oo 

g = f (id-py 1 = J2 pn > 

where P n denotes the n-th iterate of the routing matrix P = (pij , i, j G { 1, . . . , d}) , 
and the series converges because, under our assumptions, the spectral radius of the 
routing matrix P is strictly less than unity. We moreover introduce an auxiliary 
Markov chain (£„) on {0, . . . ,d}, with an absorbing state and transition proba- 
bilities pij for i 6 {1, . . . , d} and j G {0, . . . , d}. For j e {1, . . . , d}, we consider 
Tj = inf{n > : £ n = j} with the convention that inf = +oo, and we denote 

Qij = Pifa < +oo) for i, j G {1, . . . , d}. 
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2.1. Explicit Lyapunov functions. For 7 = (71,..., jd) 6 Rl, we introduce d 
vectors 74 = (7?-, . . . , 7^), i — l,...,d, with 

(2.1) H = Iog(l + Qji7i) fori,jG{l,...,d}. 

r denotes the set of all vectors 7 <E for which the following condition is satisfied: 

Definition 1. 7 £ T z/ and on/?/ if for any i — 1, . . . , d and /or any non-zero vector 
v = (v 1 , . . . ,v d ) <E M.+ with v i = 0, 

(2.2) 71 • u < sup 7^ • w 

i<j<d 

Here and throughout, u • v denotes for u, v £ R d the usual scalar product in R d . 
Our first general preliminary result is the following statement. 

Theorem 1. Under the hypothesis (A), for any 7 £ r, the function /i 7 : — > M.+ , 
defined by 

d 

(2.3) hy(x) = ^exp(^-x), ieZf, 

i=i 

satisfies the equality 

(2.4) lim sup £/i 7 (x) //i 7 (x) = — min — ^- ( — — vA. 

\x\^oo l<i<d Ga VI + 7,: / 

The proof of this result is given in Section [5j 
Remark that 

t»(tx--^)>° 

V 1 + 7i / 

if and only if 

< 7 , < £ - 1. 

If 7 G r and the last inequalities are satisfied for all i = 1, . . . , d, then the right 
hand side of (|2.4[) is negative and consequently, /i 7 is a multiplicative Lyapunov 
function for (Z(t j). 

In Section |3l we provide an example of a Jackson network with a completely 
symmetrical routing matrix, where the set T has a simple explicit representation. 
Unfortunately, in general, the explicit description of the set T is a difficult problem 
and it is of interest to give another equivalent representation of T. This is a subject 
of our next result. Here and throughout, Mi denotes the set of probability measures 
on {l,...,d} : 

Mi = {9={9 1 ,--- ,9 d ) eR d + : \\6\\i = 1} 

where \\0\\i = {O 1 ] + ■ ■ ■ + \9 d \ is the usual L 1 norm in M. d . For two vertors a = 
(a 1 ,..., a d ) and b — (b 1 , . . . , b d ) in R d we write a < b if a k < b k for all k = 1, . . . , d. 

Proposition 2.1. 1) A vector 7 = (71, ... ,7^) <E Mi. belongs to the set T if and 
only if for any i = 1, . . . , d, there exists 0i = (8j , . . . , 9 d ) £ Mi satisfying 

d 

(2-5) 7^ < 

where the vectors % = (7^, 1 < fc < d) /or j = 1, . . . , d, are defined by (|2.1[) . 
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2) Moreover, if 7, > for all i = 1, . . . , d, then 7 = (71, . . . , 7^) £ V if for any 
i = 1, . . . , d, there exists 9i = {0\ , . . . , Of) £ Mi satisfying the following condition 
d 

(2.6) 7* < y^Ojlj, whenever k £ {1, • • • , d} \ {i} and Q kl >0 

j=i 

From the above proposition it follows that the set T is open in R2_. 
The proof of this proposition is given in Section [5] 

Our following result proves that the set T is nonempty and provides an explicit 
form for some of the vectors 7 £ T. Recall that the spectral radius 1Z of the routing 
matrix P is defined by : 

1Z = inf max {pP)i/ p% 

p>0 i=l,...,d 

where the infimum is taken over all p = (p\,...,pd) with positive components 
pi > 0, . . . , pa > 0. If the matrix P is irreducible, TZ is the Perron-Frobenius 
eigenvalue and the last infimum is achieved for the left hand side Perron-Frobenius 
eigenvector p of P (see Seneta [5]). Under the hypothesis (A), the spectral radius 
7Z is strictly less than unity and consequently, the set of vectors p — (pi, . . . , pa) 
satisfying the inequalities 

< (pP)i < p, Vi£{l,...,d} 

is nonempty. Remark that these inequalities are equivalent to 

(2.7) (pP)i < Pl Vie{i,...,rf} 

Indeed, a vector p = (pi, . . . , pd) satisfies the inequalities (|2.7j) if and only if the 
vector j3 = (fli, . . . , (3d) = P~pP bas positive components ft > for all i = 1, . . . , d. 
Now since the equality /3 = p — pP is equivalent to p = j3G, then (|2.7p implies 
that < pi and < (pP)i for all i = 1, . . . , d. 

For a vector p = (pi, . . . , pd) satisfying the inequalities (I2.7|) . we define 

TZ(p) d = max (pP) i /p l 

i—l,...,d 

and we let 

(2.8) x p d = sup {x > : log(l + x) > K{p) x} . 

Theorem 2. Suppose that conditions (A) and (B) are satisfied and let a vector 
p = (pi, . . . , pd) satisfy (|2.7p . Then for e > 0, the vector 7 = (71, ... , 7^) defined 
by 

(2.9) ji = £G u /pi, for all i = l,...,d, 
belongs to the set T whenever 

n • Pi 

< e < mm — — x p . 

\<i<d G u 

Theorem Q] and Theorem [2] provide a class of explicit Lyapunov functions for 
Jackson networks. Indeed, for 7, = eGuj pi, 
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if and only if 



< e < £-(^-1 



Hence, using Theorem Q] together with Theorem [2] and the equality Gj% — QjiGu, 
one gets 

Corollary 2.1. Suppose that the conditions (A) and (B) are satisfied and let a 
vector p = (pxj • • • , pd) satisfy the inequalities (|2.7p . Then for e > 0, the function 

d d 

(2.10) h e , p {x) = h e , p (x\...,x d ) = J2 H^+eGji/pif 
satisfies 

(2.11) lim sup Ch e p (x)/h e „ (x) = — e min ( — % — -) < 

|x|->-oo ' \<%<d \Pi+e&ii pi J 

whenever the vector (eGu/pi, 1 < i < d) belongs to T and 

n / ^ • Pi ( & i 
< e < mm — — 1 

i<i<d G u \ Vi 



or sufficiently, whenever 

(2.12) < e < mm { min <j x p , £~ ( ^ " 1 



r<i<d [ [ p ' Gm yi/j 

In the above results, one can replace the vector p satisfying the inequalities (|2.7p 
by a vector j3G with /3 = (/?i, . . . , /3d) having positive components $ > 0, since 
as previously mentioned, /3 = p — pP is equivalent to p = /3G. Moreover, by 
changing if necessary e, one can assume that such a vector fj = (/3i, . . . , (3d) dehnes 
a probability measure on the set {1, . . . , d}. Then for any i — 1, . . . , d, 

d d 
Pi/Gu — S ' PjGjj/Gg = y ' PjQji — QfSi 

is the probability that a Markov chain on {1, . . . , <£} with transition matrix P and 
initial distribution /3 ever hits the state i. 



2.2. Estimates of the essential spectral radius. Now, we get some explicit 
estimates for the essential spectral radius r* . The following lower bound is obtained 
by using the large deviation results of the papers [H [6] . 

Theorem 3. Under the hypotheses (A) and (B), 

(2.13) - min i(VMl - V^l? < log K 

l<l<d (_rjj 

The proof of this Theorem is given in Section [7] 

To get an upper bound for r* e we use Theorem Q] and Theorem [2l Recall that 
under assumptions (A) and (B), the quantity 9* = — logr| is equal to the supremum 
of all 8 > for which there exists a finite set E C U\_ and a multiplicative Lyapunov 
function / : — > R+ satisfying the inequality (jl.6j) and (jl.7j) . The following 
statement is therefore a straightforward consequence of Theorem [T] 
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Corollary 2.2. Under the hypotheses (A) and (B), 

(2.14) logr,* < - sup mm (J^- v \ 

Recall moreover that Y C Rl and remark that the function 

li ( m \ 

7 — > mm Vi 

i<i<d Gu VI + 7j / 

is continuous in R 1 ^. Hence, in the right hand side of (|2.14|) . one can replace the 
supremum over the set V by the supremum over the closure Y of the set Y in Mr. 



Now the question arises of a possible equality in (12.131) . This equality holds in 
particular if the upper bound given by (12. 14)) coincides with the lower bound in 
(|2.13l) . In this respect, remark that for every i — l,...,d the maximum of the 
function 

li -> li ( r— ! Vi] 

V 1 + a ) 

over 7^ 6 R + is achieved at the point 7* = y/jii/ui — 1 and equals (^/Jm — \fv~i) 1 '■ 
Hence, if 7* = (7*, . . . ,7^) £ Y, then one gets equality in (|2.13j) . More generally, 
denote by A,; the set of all 7 € R+ satisfying the inequality 

Under our assumptions, Aj is a closed interval such that 7* £ Aj c]0, /ij/Vj — 1[ 
and clearly, Aj = {7*} for all those i = 1, . . . , d for which 

Hence, using the estimates (|2.13l) and (|2 . 14[) one will get the equality in (|2.13[) if 
there exists 7 = (71, . . . , jd) € Y with 7$ 6 Aj for all z = 1, . . . , d. In Section [3l 
we give several examples where these arguments allow to get the equality in (|2.13l) . 
Unfortunately, in the general case, the right hand side of (|2.14p is not necessarily 
equal to the left hand side of (|2 . 1 3[) (see Proposition 13.61 in Section [3] below) . In 
the general case, using Corollary 12.11 we obtain 

Corollary 2.3. Under the hypotheses (A) and (B), 

1 * ^ . f fH "i\ . n 
k>gr e < — sup e mm — < 

p,e l<i<d \pi +eGu pi J 

where the supremum sup e _ is taken over all e > and p = (pi, . . . , pa) satisfying 
(|2~7|) and (pTl"2"j) . 



3. Examples 

In this section, we give some examples for which the above results can be applied 
and in particular, equality in f|2.13[) is obtained by using Corollarv l2.2l 
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3.1. Jackson network with a branching routing matrix P. We will say that 
a matrix A — (ajj, i, j = 1, . . . , d) has a branching structure if for any i G {1, . . . , d} 
the set {j G {1, . . . , d} : aji > 0}, contains at most one element. 

Recall that under our assumptions, pa — for all i G {1, ...,d}. Hence, for 
d — 2, any routing matrix P = {pij, i,j = 1, . . . , d) has a branching structure. For 
d > 2, an example of a graph corresponding to a branching routing matrix P, with 
vertices {1, . . . , d} and ordered edges (i — > j) for i, j G {1, . . . , d} such that > 0, 
is given in Figure 1. 




Figure 1. 



Proposition 3.1. Suppose that the conditions (A) and (B) are satisfied and let 
the routing matrix P have a branching structure. Then the vector 7 = (71, . . . , 7^) 
defined by (|2.9I) belongs to the set V for any e > and any vector p — (pi, . . . , p^) 
satisfying the inequalities (|2.7p . 

Proof. We get this statement as a consequence of the second assertion of Propo- 
sition [5?T] Indeed, let a vector p satisfy the inequalities (|2.7I) . Consider a vector 
7 = (71, . . . ,7d) defined by (|2 .Q[) with some given e > 0. Then obviously, 7, > 
for all i G {1, . . . ,d}. Let us show that under the hypotheses of our proposition, 
(|2.6|) holds for any i G {1, ...,d}. If i G {1, . . . , d} is such that pji = for all 
j G {1, . . . , d} \ {i} then also Qji — for all j 7^ i and consequently (|2.6j) is trivial. 

Suppose now that for i G {1, . . . , d}, there is j G {1, . . . , d} such that pji > 0. 
Then under the hypotheses of our proposition, such an index j is unique, j =/= i, 
and 

Gki = (■!../',. V k G {1, . . . ,d} \ {i}. 
Moreover, from (|2.7[) it follows that pi > {pP)i — PjPji and consequently, 

7 f = log(l + eG ki /pi) = log(l + eGkjPji/pi) < log(l + sG^/pj) = 7* 

for all those k G {1, . . . , d} \ {1} for which Gjy > or equivalently Qm > 0. The last 
relations show that (|2.6p holds with a unit vector 6>; = (6**, . . . ,6f) where Q\ = 1 
for k j and 6^ = 1 . Using therefore the second assertion of Proposition 12.11 we 
conclude that 7 G T. □ 

When combined with Theorem [I] the above proposition implies the following 
particular version of Corollary 12.11 
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Corollary 3.1. Suppose that the conditions (A) and (B) are satisfied and let the 
routing matrix have a branching structure. Suppose moreover that a vector p = 
(pi, . . . , pa) satisfies the inequalities (12. 7|) . Then the function h e>p defined by (|2.10|) 
satisfies the inequality (|2.1ip for any 

Pi 1 1 1 
< e < mm — - 

l<i<d Gj> 



V, 



From the last statement we obtain 



Proposition 3.2. Suppose that the conditions (A) and (B) are satisfied and let 
either the routing matrix P or its transposed matrix P have a branching structure. 
Then 

1 2 
(3.1) logr* = - min — - ^/v~) 

l<i<d Ijra 

Proof. Suppose first that the routing matrix has a branching structure. Then by 
Corollary 13. 11 

Pi Vi 



(3.2) logr* < —£ min 

l<i<d\Pi+£Gii p. L 

for any e > and any vector p — (pi, . . . , pd) satisfying the inequalities (|2.7|) . 
or equivalently (see the remark below Corollary 12. ip for any p = /3G with j3 = 
(/3i, . . . , Pd) having strictly positive components /3j > for all i = 1, . . . ,d. By 
taking the limits as — > for some indices i £ {1, . . . , d} one gets (|3.2p also for 
any vector p = fJG with ft > and pi > for alH = 1, . . . , d. Using again the 
equivalence between p — f3G and j3 = p — pP 7 these arguments prove the upper 
bound (|3.2p for any vector p = (pi, . . . , pd) satisfying the inequalities 

(3.3) < p, and [pP) l < p h Vi=l,...,d. 
Remark now that according to the definition of the traffic equations (jl.2[) . 

Vi = (vP)i + \i > (vP) it Vi = l,...,d. 

If a vector p = (pi, . . . , pd) satisfies the inequalities (|3.3p . then for p — {pi, . . . , pd) 
with pi — JviPi, by Schwarz inequality, 



(pP)i = ^oPi Pii - \ u2 v oPoi , PiPi* ^ \AW = Pi! = 1, . . . , d, 



and consequently one can replace the quantities pi at the right hand side of ()3.2p by 
Pi = \Jvipi, (recall that Vi > for alH = 1, . . . , d, hence pi > for alH = 1, . . . , d). 
The resulting inequality 



log r* < — £ min 



Mi 



l<i<d \ yJViPi +£G U 

with 

\/Pi / r— I \ 

£ = ~r~ WWi-y/vi) 

l<i<d \jh 

provides the following upper bound 

^/07 . 
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Moreover, if a routing matrix P has a branching structure, then for any i € 
{1, . . . , d}, either pji = for all j € {1, . . . , d} and consequently, 

d 

?j GjjPji =0 < Gu, 

3 = 1 

or else there is a unique j € {1, . . . , d} such that Pji > and consequently, 

d 

^ ] GkkPki = GjjPji — Gji < Gu. 
k=l 

These relations show that the vector p — (pi , . . . , Pd) with pi — Gu satisfies the 
inequalities (|3.3j) . Using (|3.4|) with this vector p we obtain 

1 2 
logr* < - mm — (^/Jm - y^) . 

±<i<d (ju 

The last inequality combined with (|2.13l) proves (13. ip . 

Suppose now that the transposed matrix *P has a branching structure, and 
let us show that in this case, the equality (|3.1[) also holds. For this we apply a 
time reversing argument to the Markov process (Z(t)). The time reversed Markov 
process {Z(t)) is generated by 

my) - E %v> - /(»)). y e z ? . 

with 

q(y,z) = n(z)q(z,y)/x(y). 

A straightforward calculation shows that this is also a Jackson network but with 
different parameters: the arrivals at the i-th queue are Poisson with parameter 
= ViPio , the services delivered by the server are exponentially distributed with the 
same parameter pi — pi as for the original Jackson network (Z(t)), and the routing 
matrix (pa,i,j = 0, ... ,d) is given by p i0 = A, /', ;md /),, = Vjp^/vi for i,j e 
{1, . . . , ii}. Under our assumptions, the time reversed Markov process (Z(t)) also 
satisfies the conditions (A) and (B) with the same solution i = 1, . . . , d) of the 
traffic equations and the same stationary probabilities (tt(x); x e Z^). Moreover, 
for any finite subset E c , letting 

re = inf{i > : Z{t) E E} and t e = mi{t > : Z(t) e E} 

one gets 

V x (Z{t) =y,T E >t) = 7r(y)¥ v (Z(t) = x, t e > t)/n{x), \fx,yeZ d + \E 

and consequently, the essential spectral radius of the time reversed Markov process 
{Z(t)) is the same as for the original Markov process (Z(t)). If the transposed 
matrix P has a branching structure, then the routing matrix P = (pij, i,j = 
1, . . . , d) has the same property and consequently, the above arguments applied to 
the time reversed Markov process (Z(t)) prove the equality (|3.ip . □ 
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3.2. Jackson networks with a completely symmetrical routing matrix P. 

Now we consider a Jackson network having a completely symmetrical routing matrix 
P = (pij, i,j = 1, . . . , d) with pij ■ — p < l/(d — 1) for all i ^ j, i,j G {1, . . . , d}. 

Then Q LJ = p/(l - (d - 2)p) = q for all i ^ j, i, j G {1, . . ., d}, where < q < 1. 
The following proposition provides an explicit form for the set T in this case. To 
formulate this result, it is convenient to introduce the function 

y( x = >p maxKKrf log(l + q-fj) - log(l + qjj) 

l7l '"" 7dj " ^ log(l + 7i)-log(l + g7i) 

for 7 G satisfying jj > for all j = 1, . . . , d (note that for such a 7, since 
q < 1, then log(l +7,-) > log(l + qjj) for all j = 1, . . . , d and the above quantity is 
well-defined) . 

Proposition 3.3. Suppose the conditions (A) and (B) are satisfied and let p^ — 
p < l/(d — 1) for all i ^ j , i,j G {1, . . . , d}. Then 7 = (71, . . . , 7^) G T if and only 
if 7i > for all i G {1, . . . , d} and 

(3.5) S( 7l ,...,7 d ) < 1. 

Proof. For any 7 G any i G {1, . . . , d} and a non-zero vector v = (v 1 , . . . , v d ) G 
with v l = 0, letting \v\ = vj > 0, the inequality (|2.2p becomes 

\v\ log(l + <37i) < max (w- 7 log(l + 7-,) + (\v \ - v r )\og(l + qjj) ) 

l<j<d\ / 

or equivalently, 

(3.6) \v\ log 1 + g7 ' < ^'log 1 + 7i for some j G {1, . . . ,d\ 

1 + qjj 1 + g7i 

Since q < 1, the inequality (|3.6[) is trivially satisfied when 

7j < max7j 

Thus, 7 G T if and only if (|3.6[> holds for any i G {1, . . . , d} such that 

(3.7) 7, = maxTj 

3 

and for any non-zero vector v G with i/ = 0. 

Consider now a vector 7 = (71, ... , 7^) with ji > for all i G {1, . . . , d}. If 
7 G" r, then using the above arguments it follows that for some index i G {1, . . . , d} 
satisfying the equality (|3.7p there is a non-zero vector v G Ml with v % = such that 

Mlogi±92i > v *iogi±T*. foraU/cG 
1 + ?7fc 1 + 97fe 



and consequently, 

maxj log(l + QTj-) - log(l + 97^) 



|u| > v k for all fc G {!,..., d}. 



log(l + 7fc)- 1 °g( 1 + 97fe) 
Summing these inequalities proves that for such a vector 7, (|3.5p fails to hold. 

Conversely, suppose that 7 G T. Then (|3.6p holds for any index i G {l,...,d} 
satisfying the equality p. 71) and for any non-zero vector v G with = 0. 
From (I3.6P it follows that 7 is non-zero. Moreover, let i G {1, . . .,d} satisfy (|3.7p . 
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Then for any k £ {1, . . . , d} \ {i}, using the inequality Q3.6j) with a unit vector 
v = (v 1 , . . . , v d ) such that v k = 1 and w J = for j ^ k, one gets 

gmax7 3 = q-f t < 7 fc , 
j 

and consequently, 7fc > for all k £ {1, The quantity £(71, . . . , 7d) is 

therefore well-defined and equal to for v = (v 1 , . . . , v d ) £ R.1 given by 

j maxi<^< d log(l + (77,) - log(l + qjj) 

V ^ log(l + 7i )-log(l + g7i) ' 1 "'" *' 

If |u| = £(71, . . . ,7d) = 0, then (|3.5p obviously holds. Otherwise, using again p.6[) 
with such a vector u and with any i satisfying (|3.7[) gives 

2j(7i, 7d) max log < max log lor some j £ {1, . . . , a\ 

\<i<d 1 + qjj l<i<d 1 -(- q^j 



which proves p.5[) . 

□ 

Remark that for a completely symmetrical routing matrix P, 

(rf-l)p 2 



Hence, when combined with Theorem[TJ the above proposition implies the following 
statement, similar to Corollary 2.1. 



Corollary 3.2. Under the hypotheses of Provosition HOI the {unction 

d d 

h y (x) =^exp(^-a;) = ^(1 +7i) xi (l + ?7i) E ^ 

i=l i=l 

satisfies 

Ch 7 (x) ( {d-l)p 2 \ . ( Mi \ . n 

lim sup — — — — = — 1 — — mm 7, v: < 

[xKo ; /^(x) V l-(d-2) Py / i<<<d Vl+ 7l V 

whenever (|3.5|) and < 7$ < — — 1 for all i = 1, . . . , d. 

Note that (|3.5[) is satisfied for any vector 7 G K. d such that 71 = . . . = 7d > 0, 
so that the set of vectors 7 £ M. d satisfying both (|3.5|) and < 7, < — — 1 for all 

Vi 

i = 1, . . . , d is nonempty. Using therefore Theorem [3] and Corollary 12 . 21 we obtain 
Corollary 3.3. Under the hypotheses of Provosition [X51 

jd-iy \ . ^_ 2 , 
1 ' i-(d-2) P J ^(VR-^ < ^ 

< J^g W min 7i (-*--„/| < 

~ Ki<d VI 



1— (d— 2)p/ 7 l<i<d VI + 7 

where the supremum is taken over all 7 g T, or equivalently, over all 7 = (71, . . . , jg) 
with 7j > /or all i = 1, . . . , d such that inequality (13. 5[) holds. 
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Thus, if the conditions of Proposition 13.31 are satisfied and 



(3.8) 
then 
(3.9) 



sup min 7; ( — 

7£ r i<i<d VI 



Mi 



logr* 



7i 



min (JJTi 



min (y/JIi- 



1 - (d— 2)p J \<i<d 

that is, relation (|3.1I) again holds. The following statement gives some simple 
sufficient conditions for the equalities (|3.8I) and (13.91) 

Corollary 3.4. Suppose that for some Iq G {1, . . . , d}, 



(3.10) 
and 
(3.11) 



min (yfjTi- y^i) 



'Mi 



mm 

Ki<d 



Mi 



/Mi 



/Mio V" l ° - 

TTierj under the hypotheses of Proposition \3.3[ Q3.9P holds. 

In particular, (|3.9p ZioZds z/ one 0/ i/ie following conditions is satisfied : 

(i) /J,i/vi= fij/vj for all i, j G {1, . . . , d}, 

(ii) t/iere is io G {1, . . . , d} such that [ii > /ii and 1^ < Vi a for alii G {1, . . . , d}. 

Proof. Here, as noted above, any vector 7 = (71, . . . , 7^) with 71 = . . . = 7^ > 
belongs to the set T. Hence, by Corollary 13.31 the equality (|3.9I) holds if 

Mi 



- v, 



(3.12) sup min t (V 

t>o i<*<d \l + t 

Recall that the maximum of the function t G 



= min {y/JTi- v^) 

l<i<d 



M i 



1 + t 



is achieved at 



the point 7* — \J [lijvi — 1 and equals (y^/II 
(f3~T2"f holds if and only if 



1 >ia 



M^o 



Hence, assuming (|3.10[) . then 

Vi€{l,...,d}. 



Since 



* / Mi \ / M» ^i A / , , — \ 

the last inequalities are equivalent to (|3 . 1 1 [) . 

Now if condition (i) is satisfied, consider ig such that mini<j<d z/, 
(|3.10p is satisfied. Using ^ = iii Q Vilv iQ , we get 



then 



mm 

Ki<d 



Mi 



'Mio 



' I 1 — 
= mm — {y[Wi< 



Ki<d V. 



so that p. lip holds, hence also (|3.9p from the first part of the proof. 

Finally, if condition (ii) is satisfied, then io clearly satisfies (|3.10l) and (13. 1 1[) . so 
that (|3.9p again follows from the first part of the corollary. □ 



Remark that (ii) is in particular satisfied if [ii — [ij for all i, j G {1, . . . , d}, or if 
Vi = Vj for all i, j G {1, . . . , d}. 
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Our following result is a necessary and sufficient condition for the equality 13.8|) 
Denote 



m = min {J~\h - Jul 

0<i<d 



,2 



and consider for i G { 1 , . . . , d}, 



A, = <! t e M.+ : t { — Ui } > m 



A + t 

A straightforward calculation shows that A, = [a^ &,] with 

_ Hi-Vi-m- J + Vi - m) 2 - Aviin 



a i - o 



and 



Moreover, 



_ Mi - V% - m + VifJ-i + V_i - m ) 2 - jVjjli 

2^ 



VTkNi -ieA,c{fel+: t (^- t - Vi) > 0} 



0.^-1 



and consequently 

h > y/in/vi - 1 > ai > 0, 
where &i = a, = y \i%jvi — 1 if and only if (a/mI — \f~v~i) 2 — m. We put 



a = max a.; and 7.; = min{fa;.a| 

Ki<d 



for i — 1, . . . , d. 



Proposition 3.4. Suppose that the conditions of Proposition \3.3\ are satisfied. 
Then (|3.8I) holds if and only if £(71, . . . , 7<j) < 1- 

Proof. Indeed, suppose first that (I3.8|) holds and remark that for any i 6 {1, . . . , d}, 
since J/i > 0, the function i (7^— va ) — > — 00 as t — > +00. These functions being 



continuous on K+ , it follows that the function 



mm 7, 



Mi 



l<i<d V 1 + 7,- 

attains its maximum over the closure T of the set T at some point 7 € I\ Moreover, 
relation (13.81) proves that 

M* A • / / — / — \2 de f m 



min 7j ( — 

l<i<d Vl + 7 i J 0<i<d 

from which it follows that 7^ €E A^ and consequently, 7^ > for all i € { 1, . . . , d}. 
The quantity £(71, . . . ,7d) is therefore well defined and by Proposition ^. 3[ 

^(li, ■ ■ ■ ,ld) < 1- 

To prove that £(71 , . . . , 7^) < 1 it is now sufficient to show that (71 , . . . , 7^) achieves 
the minimum of the function £(71, . . . , jd) over (71, . . . , 7^) G Ai x • • • x A<j. For 
this let us notice that this function is continuous on the compact set Ai x ■ • • x A<j 
and hence attains its minimum on this set at some point 7* = (7* , . . . , 7^). 
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If S(7i, ... ,7^) = 0, then from the definition of the function £ it follows that 
7* = maxj 7* for all j = 1, . . . , d, that is, the intervals Aj, i = 1, • • • d, have some 
common point t = 7* = • • • = 7^ . But in this case, 

7 ^ . ^ def ~ 

mm bi > t > max cii = a 

i l<i<d 

from which, using the definition of the vector 7, it follows that 71 = . . . = 7^ = a 
and consequently, also £(71, . . . , 7^) = 0. 

Suppose now that £(71, • • ■ , 7*i) > and let us show that in this case, 7* = 7. 
Indeed, in this case, from the definition of the function £(71, . . . , jd) it follows that 
7* < max., 7* for some j = 1, . . . , d. Moreover, 

(3.13) max 7* = max at d = a. 

l<i<d l<i<d 

because otherwise, one could find some e > for which the vector 7' = (7J , . . . , 7^) 
given by 

, _ \lj- e if 7* = maxi< l < (i 7*, 

3 if lj < maxi<,<d 7* 

belongs to the set Ai x • ■ ■ x and satisfies £(71, . . . , Y d ) < £(71, ■ • ■ , "id)- Remark 
now that the following two assertions are equivalent : 

(i) (71 , . . . , 7 d ) £ Ai x • • • x A d and maxj 7, = a 

(ii) dj < 7j < min{bi,a} = 7} for all j G {1, . . . , d}. 
Moreover, for any point 7 = (71, . . . , 7^) satisfying the inequalities (ii), 

A log(l + qa) - log(l + qjj) 



£(7i,--- ) 7d) 



log(l + 7j) - lo g(! + Tii) 



When combined with (|3 . 1 3p . these remarks show that the point 7* = (7^, ... ,7^) 
achieves the minimum of the function 

log(l + qa) - log(l + g 7j -) 
^ log(l+7 J )-log(l + g 7j -) 

over the set [a 1,71] x • • • x [0^,7^]- The function 

^ log(l + go) ~ log(l + qt) 
log(l + t) - log(l + qt) 

being decreasing on ]0,a], from this it follows that 7* = jj for all j <E {1, . . . ,d} 
and consequently, 

£(7i,...,7d) = £(71, ■ • . ,7d) < £(71,..., 7d) < 1. 

Conversely, suppose that £(71, . . . ,j d ) < 1 an d let us prove the equality 
We know from Section 2.2 that 

Mi 



min (yTH - V^) 2 = sup min 7* (- 

l<i<a ^,cTl?d Ki<d \ 1 



7 gR» ^--^^ >J - 1 7i 
( Mi 

> sup mm 7i 
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Moreover, since % £ A* for alii £ {!,...,<£}, 



min % ( * - k i ) = min (y^ZJ - y^)^ 



l<i<d V 1 + 7j / l<i<ti 

To get p. 81) it is therefore sufficient to show that 7 € I\ If £(71, . . . ,7<j) < 1, then 
7 € r by Proposition 13.31 Suppose now that £(71, . . . ,7d) = 1. Then clearly 

min 7^ < max 7$, 

l<i<d l<i<d 

and letting 

/ £ \ |7i _£ ir 7« = maxi<i< d 7i, 
1 7i otherwise, 

one gets 7i(e) > for all i £ {1, ...,d} and E (71(e), . . . ,7^(2)) < 1 for all e > 
small enough. By Proposition l3.3[ it follows that (71(e), . . . , 7d(e)) £ T for all e > 
small enough and consequently, letting e — > we conclude that 7 6 I\ □ 



The last result of this section provides an example where (|3 . 8[) fails to hold. This 
example shows that unfortunately, in general, the left hand side of (|2.13j) and the 
right hand side of (|2.14[) are not necessarily equal. 

Proposition 3.5. Suppose that the following conditions are satisfied : 

(i) d > 3; 

(ii) Ai + . . . + A<j = 1 and = Ai < A, for all i £ {2, ... , d}; 
(hi) y/JTi — y/vl = t > for all i £ {1, . . . , d}. 

Then under the hypotheses of Proposition \3.3l for any p > small enough, there is 
t p > such that for t > t p , the inequality (I3.8[) fails to hold. 

Proof. By Proposition 13. 4[ it is sufficient to show that 

(3.14) lim lim > 1. 

p— >0 t— >oo 



Remark that under the hypotheses of Proposition 13.51 a, = hi = yfMj vi — 1 and 
consequently, 

it = vV»M - 1 = 



for all i = l,...,d. Moreover, a straightforward calculation shows that for any 
i = l,...,d, 



p + 1 \ 1 + p — dp ^— ' / p + 1 \ 1 + p — dp 



Since under the hypotheses of our proposition, Ai = < Ai for all i £ {2, . . . , d}, 
the above relations show that max.; % = 71. Using the definition of £(71, ■ ■ ■ , 7d) 
we conclude therefore that 

yr -x r y log ((1 + qt/^)/(l + gt/y/vt)) _ ' log (vj/vi) 
t nn l7 i,..., 7d j log((l+t/v^)/(l + 9t/v^)) &21og(l/g) 

where q = p/(l + 2p — tip) and 

Vijv\ = 1 + (1 +p— dp)Xi/p, Vi = 2,...,d. 
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Hence. 



lim E(7i,...,7 d ) 

t— >oo 



and since Xi > for i — 2, • • • , d, 



21og(| 



lim lim £(71, . . . , 7 d ) 

p— >0 t— too 



2-d 



d-l 



Under the hypothesis (i), the last relation proves (|3 . 14|) . 



□ 



3.3. Jackson network with three nodes on a circle. Consider a Jackson net- 
work with three nodes (d = 3) and a routing matrix 

/0 p q\ 

(3.15) P=\q p\ with < p < q < 1 such that p + q < 1. 
V 9 0/ 

(see Figure 2). Here, as a consequence of Corollary 12.21 and Theorem [3] we get 

P 




Figure 2. 



Proposition 3.6. Suppose that a Jackson network with three nodes and a routing 
matrix (13. 15)) satisfies conditions (A) and (B). Then 

(3.16) - min (y/Ji - y^) < 

, * , l-p 3 -q 3 -3pq . / /Xj 

logr e < — sup mm rl 



\ — pq t>o i<i<d \l + t 

If moreover the equalities (|3. 10[) and (13. lip /io/c! /or some io G {1, . . . , d}, then 

(3.17) logr* = min - v^i) • 

1 — pq % 

In particular (|3.17p holds if at least one of the conditions (i) or (ii) of Corollary 
is satisfied. 
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Proof. Indeed, a straightforward calculation shows that 



I Q +P P +q 

G = f (Id-P)- 1 = - = ___ U 2 + g 1-pq q 2 +p 



1 — p 3 — q 3 — 3pq 



q 2 + p p 2 + q 1 — pqj 



The first inequality of f|3 . 16[) is therefore a straightforward consequence of Theo- 
rem [3l By Corollary 12.21 to prove the second inequality of (|3.16p it is sufficient to 
show that for any t > 0, the vector 7 = (t,t,t) belongs to the set T. For this let 
us first notice that under the hypotheses of our proposition, the matrix of hitting 
probabilities Q = (Qij, i,j = 1, 2, 3) is given by 

. (l—pq q 2 + p p 2 + q^ 

Q = - — — ( p 2 + q i — PQ q 2 +p 

K q 2 +p p 2 + q l-pqj 



1 — pq 



Without any restriction of generality we can assume that p < q. Then 

Q31 = Q12 = Q23 < Q21 = Q32 = Q13 < 1 
and consequently, for 7 = (t,t,t) with t > and any v = (1)1,1)2,1)3) € Mr, with 
Vi = and (^1,^2) 7^ (0, 0), one gets 

71 -d = v 2 log(l + Q 2 it) + v 3 log(l + Q 31 t) 

v 2 log(l + t) + v a log(l + Q 32 t) = 72 • v if v 2 > 0, 



< 

from which it follows that 



t>3 k>g(l + t) = 73 • V if W3 > W2 = 0. 



71 ■ u < max 7^ • v. 

3 



Permuting indices shows that for any i £ {1,2,3} and any non-zero vector v — 
(vi,v 2l v 3 ) e R\ with Vi = 0, 

7i • v < max jj ■ v, if 7 = (t, t, t) with t > 

i 

Hence, for any i > 0, the vector 7 = (i, t, t) belongs to the set T and consequently, 
by Corollary 12.21 the second inequality of (|3.16|) is also verified. The first part of 
our proposition is therefore proved. The second part of Proposition l3 . 61 follows from 
(|3.16[) by using the same arguments as in the proof of Corollary 13.41 □ 

4. Background 

For a given A C {1, . . . , d}, denote A c = {1, . . . , d} \ A and consider the sets 

R A,d dej ^ xeM d . x j > Qj y- g A | and 

B A = f ^a = (a 1 ,...,a d )€M d : a 1 < log(^p tf e Q ' + Pl0 ), Vi ^ A 
For (3 € R+' and i,j G {0, 1, . . . , d}, 0, we define 

m ij W) = Pije~ pZ + Yl PihPhh ■ ■ ■ PJnj ex P I ~P % ~ Y P 

n>lji,...,j n eA" \ 

The following result provides a suitable homeomorphism from the set M/^' d onto 
B\, this is a straightforward consequence of Proposition 8.1 of the paper [j. 
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Proposition 4.1. (Proposition 8.1 [4]^ Under hypothesis (A), 

- for any Ac {1, . . . , d} and € R^' d , the system of equations 

ft = a\ for i g A, 

ft =log{j2 d 3=lPlJ e a3 - aZ + Pl0 e~ a ') , for lE {l,...,d}\A 

has a unique solution a = a\((3) g Ba : 

f ^(ft = ft forieA, 
(4 - 1} | a'CS) = log(Ei e A<(^)^+<(^)) VzeA c , 

- the mapping (3 a\(f3) determines a homeomorphism from onto the 
set B\; 

- the function R(oca((3)) is strictly convex in R^;' d . 

This result will be used to investigate the different Laplace transforms of the 
jump distribution on the different "faces" of the space Z+. For Ac {1, • • ■ , d} and 
a g K d , the Laplace transform of the jump distribution corresponding to the face 
A is defined by 

d / d 

Ra{o) = f ^(e^ -l) + Y f N[T f Pike a "-^ +p j0 e-^ -1 
3=1 jeA \k=i 

and for A = {1, . . . , d}, we denote 

d d I d 

R(a) = f = E^^-^ + E^ Ep^-^+p j0 e- Q, -l). 

3=1 3 = 1 \k=l ) 

As a consequence of the above proposition one gets the following statement . 

Lemma 4.1. Under hypothesis (A), for any i g {1, • • • ,d} and s g ] — 1, +oo[ the 
system of equations 

d 

(4.2) e a% = l + s, e a) ^ft-Z+ft-o, j g {1, . . . , d] \ {i} 

k=l 

has a unique solution a — a(s) — (cr^s), • • • ,a d (s)) given by 

a j (s) = log(l + Q jiS ), j €{!,..., d}\{i}. 
Moreover, for any A C {1, ■ • • ,d}, this solution satisfies the equality 

«,(„(.)) . _£.(„__*_!,<„). 

Proof. Indeed, for /? = (ft, . . . , j3 d ) with ft — log(l + s) and /3 J = for j ^ i, one 
gets 

™£ } (/?) = Qji and mg } (/3) = 1 - Q J4 
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for all j £ {1, . . . ,d}. Hence, the first assertion of Lemma |4J] is a straightforward 
consequence of Proposition 14. II Moreover, for any A C {1, . . . , d}, 

R A (a(s)) =J2*j (^ (S) - l)+IH ]>> fe e afc(s) - Q,(s) + Pj0 e- a '^ - 1 l A (i) 

3=1 \k=l J 

d / d \ 

j=i s \fc=i / 

d / d \ 

= s^AjQjj + — ^— I ^PikQki - 1 J 1a(«) 
j=i s \fe=i / 



The last equality combined with the relations 

d , d 



y^AjQjj — — ^ AjGji — — ^- and 1 -^^PikQki — tt - 

j=l 3=1 fc=l 



proves the second assertion of Lemma 14.11 □ 

5. Proof of Theorem [TJ 

We begin the proof of this theorem with the following lemma. 

Lemma 5.1. For any 7 £ Mi and 1 < i < d, the function fi(x) = exp(7i • x) 
satisfies the equality 

(5.1) Cfi(x) = ^- (ut - l {xi>0}r a-^ f t {x), x£7L\ 

Proof. A straightforward calculation shows that for any a = (a 1 , . . . , a d ) £ R+, 
the exponential function f a (x) — exp(a • x) satisfies the equality 

Cf a (x) = RA(x)(a)fa(x), 

where for x £ Mi, we denote by A(x) the set of all j £ {1, . . . , d} for which x J > 
and for A C {l,...,d}, 

d I A \ 

R A (a) = J2 " 1) + X>i EPf e ° l " ai +P^ aj - 1 



Oj k e 

3=1 jeA \k=i 



Furthermore, by Lemma 14. 1 L from the definition of the vector 74 it follows that 
a = (a 1 , . . . , a d ) = 7I is the unique solution of the system (|4.2j) for s = 7^ and 

i?A(,)(7f) - g~ (** ~ 1{iBi>0} T+yU 
The equality (|5.ip is therefore verified. □ 

Now we are ready to complete the proof of Theorem [T] For the function /i 7 
defined by (|2.3[) , Lemma 15.11 proves that 



(5.2) Ch 1 (x)= E 77- ( ^ " -, ,\. ) exp(7l ■ x) + ^ TT 1 exp(^ • x 
- m ^ x > \ Vl ~ TI~ ) M^H E ^-7»exp(7l -x). 

i£A(x) 



22 



I. IGNATIOUK- ROBERT AND D. TIBI 



To get the inequality 

li ( V 

l sup L,n-y \x) I n-y (X) v max 

\x\ 

it is therefore sufBcient to show that 

max^ A(l ) exp(7f • x) 



(5.3) lim sup £/i 7 (x)//i 7 (x) < max -^—\Vi 

l<i<d Gii V 1 + 7i 



max 



A(s) exp(7, • x) 



as |x| — > oo, 



or equivalently that 



(5.4) lim exp max ji ■ x — max 7, • x ) =0. 

s-oo yi^A(a;) ieA(a:) / 

The last relation follows from the definition of the set T. Indeed, using the inequality 
(I2.2[) with v = x/\x\ for an arbitrary x £ R+ \ {0}, one gets 

from which it follows that 

v X ^. X ^. X 

max 7I • - — - < max 7I • - — r = max 7I • - — - 

i^A(x) \x\ , 1 ./ |x| jeA(x) \x\ 

for any non-zero x £ R+. The function 

x — > max 7I ■ x — max 7$ • x = max 7! • x — max 7} ■ x 

i^A(x) i£A(x) igA(x) j=l,...,d 

being upper semi-continuous on the compact set Sf. — {x £ : |x| = 1}, from 
the above inequality it follows that 

(5.5) max 7^ • x - max 7^ • x < -5\x\, Vx £ Rl \ {0} 

i^A(x) ieA(x) 



with some 5 > 0, and consequently, (|5.4[) holds. The inequality (|5.3|) is therefore 
proved. Moreover, (|5.2p and (|5.5p applied for x £ R+ with A(x) = {i} prove that 

limsup Ch~ f (x)/h 1 {x) = ( z/j 



A(x)={i} Gu \ 1 + 7» 

Using this relation together with (|5.3p one gets (I2.4[) . 



6. Proof of Proposition 12. II 
We begin the proof of Proposition 12.11 with the following lemma. 

Lemma 6.1. For u\,--- ,Ud £ R d . the following two properties are equivalent: 

(1) for any v £ R^ \ {0} ; there exists some i £ {1, • • • , d} such that Ui ■ v > 0, 

(2) there exists some 9 = (6>\ • • • , 9 d ) £ M\ such that Y?j=i > °- 

Proof. It is straightforward that (2) => (1), since for 9 satisfying condition (2) and 
for any v £ R^ \ {0}, 

d d 

< v ■ 9 j Uj = 9 J Uj ■ v 
so that one of the non-negative terms of the last sum needs to be positive. 
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To prove the converse, assume that (2) is not satisfied, so that for any £ A4i, 

d 

Y^° ju 3 $ ]0,+Oo[ d 
.7=1 

This means that the two convex subsets of M. d given by the open orthant ]0, +oo[ d 
on one hand, and the closed convex cone C generated by vectors U\,-" ,Ud on 
the other hand, that is, C = |X)j=i ®j u j > ® — (^Ij ' * ' i @d) € are disjoint. 

Then by Hahn-Banach theorem, there exists some hyperplane separating these two 
convex sets, that is, there exists some v € R d \ {0} and some c € K such that 

]0, +oo[ d C {x e M. d : x ■ v > c} and C C {x e R d : x ■ v < c}. 

Note that the first inclusion extends to the closed orthant [0,+oo[ d . Now since 
the zero vector is both in C and in the closed orthant, the constant c must be 
zero. The first inclusion, extended to [0, +oo[ d and applied to the canonical vectors 
ej for i = 1, • • • , d, yields that v has non- negative components. And the second 
inclusion above implies in particular that Ui ■ v < for all i, proving that (1) is not 
satisfied. □ 

We are now ready to complete the proof of Proposition 12.11 
For 7 € Ml, the condition 7 £ T, described by the inequalities (|2.2[) . says that 
for any i € {1, • • • , d}, the property (1) of the lemma is satisfied, with d— 1 in place 
of d and with, as vectors Mi's, the d — 1 projections on Ri 1 '"' > d }\W Q f the vectors 
7j — 7j, j 5^ i. The lemma thus proves that 7 € T is equivalent to existence for 
each i = 1, • • • , d, of some 0j G A4i satisfying 0| = and for all k € {1, . . . , d} \ {i}, 

d 

(6-1) 7 f < 

3=1 

It is straightforward that the condition 9\ = can be removed. The first part of 
Proposition 12 . 1 1 is therefore proved. 

Suppose now that 7 > 0, and that for any i — 1, . . . ,d, there exists some Q{ = 
(9\ , . . . , 9f) E Mi, satisfying the inequalities (|6.1I) for all those indices k for which 
Qki > 0. The inequalities (16.11) being strict, without any restriction of generality, 
one can assume that Of > for all i, j € {1, . . . , d}. Then for all i, k G {1, . . . , d} 
for which Qki = 0, 

d 

7 , fc = < 0?log(l+7*) = Ohk < 
The inequalities (|6.ip hold therefore for all i, fc € {1, which ensures that 

7 e r. 

7. Proof of Theorem [2] 

Suppose that the conditions of Theorem [5] are satisfied. For the vector 7 = 
(71, . . . , 7 d ) defined by ([2T5]) . it follows from (|2TT|) that 

7^ = log(l+ eGji/p,) > 0, Vi,j€{l,...,d}. 

Since 7^ > for all i = 1, • • • ,d, then from the second assertion of Proposition 12. 11 
7 G T if for every i, (|2.6p is satisfied with some vector 0i e M.\. Let i € {1, • • • , d} 
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and k ^ i be such that Qki > 0. Then, there is j 6 {1, . . . , d}\ {i} such that pji > 
and consequently, (pP)i > 0. Letting Q\ = PjPji/{pP)i for j = 1, . . . , d we obtain 



7i 



log(l + £ G fci / Pi ) = log l + < £ £ 



(P-P)i <sp c GkjPji = (pP)i y> gj c G kj 
Pi {pP)i Pi £i 1 Pj 



(7.1) < n(p)J20ie^-. 

Assuming now that 



r, • Pi 

v < e < mm — — x p , 

l<i<d (jn 

one gets 

< £— ^ < x„ for all J G {1, • • • , d}, 

Pj 

where the left inequality is strict at least for some j £ {1, • ■ ■ ,d} with pj t > 0, 
because Qki > implies that Gk% — Y^j=i GkjPji > 0. It then results from the 
definition of x p that 



ft(p) e — ^ < log 1 + £— for all j £ {1, ■ • ■ , d}, 
Pj \ Pj J 

where the inequality is strict at least for some j with Q\ > 0. The last inequality 
combined with (|7.1|) proves that 



3=1 V PW J=l 



The condition (|2.6I) of the Proposition 12. II is thus satisfied and therefore, 7 G T. 

8. Proof of Theorem [3] 

To prove Theorem^ we use the equality (jl.101) and the explicit representation 
of the sample path large deviation rate function JVt] (4") obtained in [H [5] . Recall 
that the family of scaled processes Z e (t) — eZ{t/e), t € [0, T] satisfies the sample 
path large deviation principle (see [TJ [2l [4l [5] ) with the good rate function 

J J T L(<j)(t), 4>{t)) dt if 4> : [0, T] -> is absolutely continuous 
+00 otherwise 



I[0,T] (<P) 



where the local rate function L(x,v) is given by the formula (see [4]) 
L(x,v) d = sup (a-v-R(a)), Vv £ R d , x E %\. 

As above, a ■ v denotes here the usual scalar product of a and v in R d , 

d d d 

R(a) = f ^^^pye^-^+pftje-^-i^+^Aife '-!) 
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for x = (x 1 ,...,x d ) G R^_, 

A(sb) = J {« G {l,...,d} : x l > 0} 
and Ba is the set of all those a = (a 1 , . . . , or) G Mr for which 

d 

e Q * < ^Py-e Qj +Pio for all i A. 
i=i 

For a constant function <j) x {t) = a; with a: G (a; 1 , . . . , x d ) G Ki, we get 

I[o,i]{<t>x) = - inf i?(a) 

tt£BA(.) 

and using (jl.lOp we obtain 

logr* > — inf I\Qi\{<t>x) = max inf R(a) 

> max inf R(a). 

l<i<d QSB {i} 

To prove Theorem [3] it is therefore sufficient to show that 

(8.1) max inf R(a) = — min — — (JlTi — \/ui) 2 - 

l<i<d aeS f ,j l<i<d Ga 

For this we first notice that for any i G { 1 , . . . , d} , 

inf R(a) = inf { R(a) a G R d , e a ' < Vp^e"' +p j0 , V j ^ i 1 

< inf jiZ(a) | a G R d , e« J = Yj*^ + PiO, V j £ * j 
Lemma 14.11 shows that the right hand side of (|8.2[) is equal to 



7« A. /'< \ _ 1 , /t- /7^2 



7*>-l &u V 1 +7iJ 

Without any restriction of generality we can assume that 

(8-3) min — (y^ - ^/Ulf = — - v 7 ^) 2 . 

l<i<d Crjj On 

To get (18.11) . it is now sufficient to show that (|8.2[) holds with the equality for i = 1. 

For a given a G 1^, it is convenient to introduce the set J(a) of all those 
j G {1, . . . , d} for which 

e = 2_^Vkj£ +p k Q- 

The proof of equality in fl8.2[) for i = 1 uses the the following lemma. 

Lemma 8.1. Suppose that the conditions (A) and (B) are satisfied and let (|8.3|) 
hold. Suppose moreover that a G £?{i} and {2, . . . , d} \ J(a) ^ 0. Then for any 
i G {2, . . . ,d}\ J(a), there exists an a £ B^iy such that J(a) U {i} C J{a) and 

f 1 2 

i?(d) < max^i?(a), - — (y^ - y^J 
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Proof. Indeed, consider the vector ji = (7/ , . . . ,7^) defined by (|2.ip with 7, = 
7* = \/ Vil v i ~ 1 > 0. Then for k ^ i, using Lemma f4.lt 

d 

(8.4) /'/,,< ' • /'„, = e 7 - , 
from which it follows that 7^ G S^j C <B{i,j}- Moreover, 

j'=i j'=i j=i 

(8.5) = l+fl- J-) 7 ? < 1 + 7* = e 7 * 

and 

(8-6) iZ(^) = £L--J±-) = 7^(VMl-v^) 2 , 

Consider now the homeomorphism j3 — > asin{0), from R^ 1 '*^ to B{i^, defined 
by (|4.1[) for A = {l,i}, and let a fir in (a) denote its inverse mapping. Then 
the equality (|8.4[) implies that a(yi) = for all k G {1, . . . , d} \ {1, i}. Suppose 
now that a G &{i} an d i G" J {a). Then according to the definition of the set #{1}, 

d 

^Pije" 3 +p i0 > e Q * 
i=i 

and a G 23{i ^} . Since the function .R(o:.n $} i s continuous, the last relation 
combined with (|8.5[) shows that for some < s < 1, the point /3 = s/3{i ) j}(7i) + 
(1 — s)/3{i.i}(o0 G '' d satisfies the equality 

(8.7) J2 PlJ e a3 i^+ Pi0 = e«U.«}^ 

j'=i 

and consequently, i G J(a.n ) i} (/?)). Moreover, j3j — for all those j G {1, . . . , 

for which (3^ ^(a) = and consequently, J (a) C J(a{ l! ji.(/3)). Finally, recall that 

by Proposition 14. 11 the function R(otJi t a(/3)) is strictly convex. Hence, 

R{a { i ti} 0)) < mm{R(a) f R($)}, 

and therefore, our lemma is verified with a — a.n,i}(/3)- □ 

Now we are ready to complete the proof of Theorem[3] By induction with respect 
to the set J (a), for any a G B^ with J(a) j& {2, . . . , d} there is a point a G B^ 
with J(a) 2 {2, . . . , d} such that 

f 1 2 

R(a) < max <i?(o:),— min — — {JJm, — y/uT) 

2<i<d (j;; 



When combined with (18.31) and (I8.6P for i = 1, the last inequality shows that 



i?(a) < max \ R(a),- — (^JH - y/vT) } = max {R(a), #(71)} ■ 
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where, as in the proof of the last lemma, 71 = (7*, . . . , 7^) is defined by (j2.1[) with 
71 = 7i = y/ '[i\jv\ — 1. Since J (a) = </(7i) = {2, ...,d} and the minimum of 
R(a) over a € M. d with J(a) = {2, . . . , d} is achieved at the point 71, using the last 
inequality we conclude that 

#(71) < R{a) < R(a). 

This proves that the minimum of R(a) over a € ^{1} is achieved at a = 71 and 
consequently, equality holds in (|8.2[) for i = 1. The proof of Theorem[3]is complete. 
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